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Yildirim et al. have performed a Monte Carlo (MC) 
study of the ferromagnetism in Ga\- x Mn x As In this 
study, based on the dilute Kondo lattice model, a realis- 
tic band structure for the host is included. We argue that 
their model is inappropriate for GaMnAs. (i) The p — d 
coupling used in their calculations is in fact too small and 
does not correspond to the realistic value of 1.2 eV.(ii) A 
realistic value leads, in contrast, to very large Curie tem- 
peratures (Tc) inconsistent with the measured ones.(iii) 
Their MC study disagrees and contradicts previous ab- 
initio calculations. (iv) The experimental data plotted in 
Fig. lb of [l[ are incorrectly used: Tc were measured with 
samples of different Mn concentration whereas their cal- 
culations performed at fixed density. 

J is often denoted N/3 in the literature, especially in 
ab-initio based studies 0, HI, Q. It is well known that 
N/3 « 1 eV in most of the II- VI materials as ZnM- 
nTe, CdMnTe 0. Its amplitude in GaMnAs was of- 
ten discussed and a value ranging from 1.2 to 3.5 eV 
(e.g. 2.5 < JS < 8.8 since Mn 2+ spin is S = 5/2) 
was proposed [ESQ- However, it is now well accepted 
that a realistic value of N/3 is close to 1.2 eV [3J, thus 
JS w 3 eV (see the expression of N/3 in [3[) should be 
used in realistic model studies. Although, Yildirim et 
al. quote this value for J, their calculations are done for 
JS = 1.2 eV since S = 1 was used. Thus, their cou- 
pling is about 3 times smaller than the realistic value. 
As a consequence, according to Fig. la of [l| the realistic 
J would lead to Tc ~ 700 K (5 times larger than mea- 
sured). On the basis that a realistic band structure was 
already taken, an agreement with experimental measure- 
ment (at x — 0.085) can be achieved within their model 
only if a non realistic small value of N (3 is assumed. For 
this value, it is not clear that they could reproduce Tc 
as a function of Mn 2+ as done within ab-initio calcu- 


lations [1, [t| [T(| • Additionally, from the ab-initio stud- 
ies (sophisticated realistic calculations) it is clearly seen 
that GaMnAs exhibits a preformed impurity band [J, [HI], 
Mn 2+ affects significantly the host band. For the small 
JS used in |l[ such effects are absent as seen in Fig. 4 
where the DOS at JS — 1.2 eV is almost unchanged 
with respect to the host DOS. Thus their J corresponds 
to the perturbative regime (RKKY). In this regime, it is 
expected that the couplings may developan RKKY tail 
that could eventually suppress the Tc [131 • Because of 
the very small clusters used in the MC simulations and 
insufficient averaging (5 configurations only), their cal- 
culated Tc may not be reliable. They have used systems 
with typically only 20 spins. As seen, for example in [Iol |. 
even with much larger systems Tc fluctuates already sig- 
nificantly from sample to sample. 

In Fig. lb of [l| a comparison between the calculated 
Tc as a function of hole density and experimental data is 
shown. By reading carefully ref. [R| from which the ex- 
perimental data were taken, one finds that the samples 
have different Mn concentration. Only the samples with 
the highest Tc ~ XAQK have x w 0.085. The comparison 
with these data cannot be done, (Tc also depends on x). 
Thus a comparison would re quir e the knowledge of the 
total Mn 2+ and hole densities [lj] . Hence, the agreement 
observed in Fig. lb of ref. l] is surprising and accidental. 
In addition as mentioned above, because their calcula- 
tions are performed in the weak coupling regime, the sta- 
bility region expected for the ferromagnetism should be 
smaller (the RKKY oscillations of the Mn-Mn exchange 
should suppress the ferromagnetic phase), in their case 
the Tc increases almost linearly up to large hole density, 
this is in contrast with the results found from ab-initio 
based studies for GaMnAs @. We conclude that their 
model is inappropriate for GaMnAs 
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